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PREFACE 


The research progress on the NASA Grant NSG 3032 is briefly reported 
here. The period covered is approximately nine months ending June 30, 1976. 
The research was initiated in September 1975. 

The report is divided into four sections as follows. 

1. Literature survey on the end wall flows in the axial 
flow compressors. 

2. Analysis of the rotor blade end wall flows in the 
compressors. 

3. 


4 . 


Experimental facility, instrumentation and technique. 
Status of the experimental work. 



I. LITERATURE SURVEY ON ENDWALL FLOWS 
IN AXIAL FLOW COMPRESSORS 

The flow near the tip region of a compressor rotor consists of secon- 
dary flows, casing boundary layer, blade boundary layer (including the 
radial migration of low momentum and energy flow) , leakage flow and scrap- 
ing vortex and hence is extremely complicated. The qualitative nature of 
end wall flows in a compressor rotor is illustrated in Fig. 1. 

A brief summary of the previous work carried out on end wall flow 
phenomena will be presented in this section with major emphasis on annulus 
wall boundary layer. 

A first serious attempt to incorporate boundary layer equations into 
compressor performance was made by Jansen (1) , who simply applied the one- 
dimensional von-Karman momentum integral equation. Although some justifi- 
cation for this procedure was provided by Stratford (2) , it is now clear 
that the application of only the von-Karman momentum integral equation is 
not sufficiently accurate, since the important effects of incremental end 
wall work and losses are missing in the crude model. Moreover, the effects 
of tip leakage as well as secondary flow were not included in this model. 

In incompressible flow through repeating stages, the wall boundary layer 
should grow to some asymptotic size and remain constant thereafter. On 
the other hand, the von-Karman equation, in the repeating stage situation, 
can only predict indefinitely growing layers until the inner and outer 
annulus layers merge. 

Recognizing the insufficiency of von-Karman momentum integral equation 
to compressor annulus wall boundary layers, Mellor and Wood (3) proposed a 
new theory for integral representation of the end wall boundary layer. 
Inherent in the derivation is a description of the role of secondary flow. 
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tip clearance flows, blade force and the wall shear stress, thus providing 
a frame work for the incorporation of each of these effects into the over- 
all problem. This theory needs empirical information such as shape 
factor, skin friction coefficient and leakage coefficient. This model is 
incorporated in performance calculations of axial compressor by Mellor 
and Balsa (4,5). Their predictions of compressor performance are good, 
except for the empirical parameters and absence of adjustable parameters. 
Hirsch (6) has applied Mellor and Wood (3) theory to predict overall per- 
formance of a single stage compressor. He has shown that the constant 
shape factor and skin friction coefficient assumptions are invalid. He 
has compared this theory with the existing data and has shown that skin 
friction coefficient takes up unusually high values at the tip of rotor 
blade. Finally, he reformulates the theory in streamline coordinates thus 
allowing for shape factor variation and the Reynolds number effects, intro- 
duces the wall skewing angle, modifies the original closure equation by 
allowing non-collateral character of the boundary layer at the blade outlet. 

Railly (7) improved the theory of Mellor and Wood (3) by deriving a 
second equation for the growth of the end wall boundary layer which relates 
the work of the mean blade force along the mean velocity with the energy 
of secondary disturbances in the boundary layer. Secondary flow theory 
is used to obtain two equations for two boundary layer parameters. Railly 
and Sharma (8) compared Railly' s (7) theory with experimental results 
downstream of a single stage rotor and have obtained good agreement between 
theory and experiment. 

Work on annulus wall boundary layers in axial turbomachines, especially 
by the Cambridge University group was reviewed in Reference 9. Recently, 
Reynolds (10) conducted experiments on axial flow research fan of Applied 
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Research Laboratory of the Pennsylvania State University. He measured 

t 

rotor exit flow up to one chord length downstream in the tip region, using 

: | 

a wedge probe and determined that the predominant flows in the end wall 
region of a compressor rotor were annulus wall boundary layer and tip leak- 
age flows. 

The blade boundary layer on a rotating blade and their interaction 
with the annulus wall boundary layer has been studied by the Pennsylvania 
State University group (11,12). Their experimental work is confined to 
axial flow inducer. The three-dimensional form of momentum integral equa- 

f 

i 

tions is derived. Eventhough this work represents one extreme case, where 
viscous effects are large, the results are applicable to axial compressors 
and other turbomachinery. 
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II. ANALYSIS OF ROTOR BLADE END WALL FLOWS 
IN COMPRESSORS 


2.1 Introduction 

The objective of the analytical investigation is to model end wall 
flows including leakage and secondary flows, casing and blade boundary 
layers and to include it in the analysis and design of axial flow com- 
pressors. So far casing and blade boundary layers were included in the 
end wall flow model. 

The approach for the analysis, the assumptions made and the important 
equations derived will be presented in the following sections. 

The analysis is a preliminary one and is being continuously updated. 


2 . 2 Nomenclature 
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Tangential direction 

Angular pitch of the annulus wall boundary layer 
region between the blade boundary layers 
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Pressure surface of the blade 
Suction surface of the blade 
Blade surface 

Radial, tangential and axial directions respectively 
Streamwise and normal directions respectively 
Edge of blade boundary layer on suction surface 
Edge of blade boundary layer on pressure surface 


Superscripts 

* Refer to the boundary layer on the pressure side of 

the blade 

Pitch average of fluctuating quantities 
' Spatial fluctuations 

" Temporal fluctuations 

2. 3 Assumptions 

The following assumptions are made in the analysis of end wall flow 


1. Flow is incompressible and steady relative to the rotor. 

2. Blade camber is small. 

3. Coles velocity profile is assumed for the blade boundary layer 
velocity in the z direction and Hager's velocity profile is 
assumed for the blade boundary layer velocity in the r direction. 

4. 6 in the blade boundary layer is assumed to be constant in the 
tangential direction. 
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2.4 Approach 

The time-averaged equations of continuity and momentum are written in 
rotating cylindrical coordinate* system (see Fig. 2). Afterwards, the 
momentum equations are integrated in 6-dlrection. The integration is 
carried in three regions: 1) blade boundary layer on pressure side, 

2) annulus wall boundary layer region between the two blade boundary layers and 

3) blade boundary layer on suction side. In integrating momentum equations 
in blade boundary layers. Cole's profile for axial velocity and Mager's 
profile for radial velocity are used. Then the momentum equations are 
integrated in radial direction within the annulus wall boundary layer. 

2. 5 Momentum Integral Formulation 

The continuity and momentum equations in cylindrical rotating coordi- 
nate system are given by, 
continuity equation: 
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Assumed velocity profiles in the blade boundary layer: 

Cole’s velocity profile is used for the blade boundary layer velocity 

in the z-direction. 

Refer to the Fig. 2, let 

0* *0-8 or 0-6. 

p s 
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Hager '8 profile is used for the blade boundary layer velocity in the 
r-dlrectlon. 
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For the present formulation, normal velocity in the blade boundary 
layer is assumed to be very small compared to the streassrise velocity in 
the blade boundary layer. 

Hence u = u cosB and u fi * u sinB. 

ZB V B 

Integration in 6-direction from blade-to-blade: 

Let q be any flow quantity, then (Fig. 2), 
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How Integrating the equations of motion in the 0-direction, and 
uiing the continuity equation in the momentum equations, we obtain the 
following equations. 

Continuity equation: 
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z-momentum equation: 
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6-momentum equation: 
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In equation (15) to (18) the starred quantities refer to boundary 
layer on pressure side of the blade. 

Let 6 w - annulus wall boundary layer thickness. 

U g ■ pitch averaged velocity at the edge of the annulus 
wall boundary layer. 

' U e ( 2 ) 

The assumption that U e is a function of z only is valid in the end 
wall region. In the other regions U g is a function of both z and r. 
Define the following annulus wall boundary layer thicknesses. 
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Multiplying continuity equation by (u z - U g ) and adding to the z- 
momentum equation and then integrating the resulting equation with respect 
to r from (r £ - 6 w ) to r^, we get the z-momentum integral equation 
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the thickness of the boundary layer on the pressure side of 


the blade. 
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III. EXPERIMENTAL FACILITY, INSTRUMENTATION 
AND TECHNIQUE 

Modifications to the compressor unit: An overall view of the com- 

pressor unit is shown in Fig. 3. A rotating traverse mechanism was 
mounted on the rotor shaft immediately after the rotor. This mechanism 
rotates at the same speed as the rotor. This mechanism along with the 
rotor was dynamically balanced at the Applied Research Laboratory of The 
Pennsylvania State University. 

The traverse mechanism is capable of giving three degrees of freedom 
to the probe used in the end wall measurements. One mode is automatic and 
the two other modes are manual. Tangential movement of the probe is auto- 
matic and is achieved by mounting the probe on a disc which is rotated at 
an angular interval of 1/10 degree by a step motor. Radial and axial 
movements of the probe are to be carried out manually. The traverse 
mechanism has six axial slots on its periphery. The circumferential length 
of these slots allows measurements in two adjacent blade passages of the 
rotor in the end wall region of the blade. The axial spacings of the slots 
give accessibility for measuring the flow. in the end wall region at various 
axial stations right from the leading edge of the rotor blade to its trail- 
ing edge. Provision has been made for measuring the hub wall and the 
blade static pressures. 

The transverse mechanism is mounted immediately after the rotor. A 
scanivalve will be used to transmit blade static pressures and hub wall 
static pressures. 

Instrumentation: The instrumentation assembly designed for trans- 

mitting the electrical and mechanical signals from a rotating reference 
frame to stationary instruments is shown in Fig. 4. This Instrumentation 


1 
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assembly enables transmission of data from a five-hole pitot probe and a 
three sensor hotwire. The assembly consists of a pressure transfer device 
(PTD 1) with a scani valve mounted in it, a six-channel brush slip ring, 
a three-channel pressure transfer device (PTD), a six-channel mercury 
slip ring and a slotted disk. The PTD 1 and PTD will be used to transmit 
rotor blade static pressures, hub wall pressures immediately after the 
rotor and pressure measurements from the five-hole probe. The mercury 
slip ring will be used to transmit signals from the hotwire. The brush 
slip ring unit will be used to drive the step motors in the scanlvalve 
and the rotating traverse mechanism. A photo cell mounted on the support 
plate of the instrumentation assembly will be used to measure the speed 
of the rotor accurately. The disk with slots will activate the photo cell 
At present the instrumentation is being assembled. The assembly will 
be connected to the rotor shaft by means of a flexible coupling to reduce 
transmission of rotor shaft eccentricity to the assembly. The instruments 
tion will be enclosed in a stationary housing to avoid flow distribution 
upstream of the rotor. 


1 
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IV. STATUS OF EXPERIMENTAL WORK 

The experimental work carried so far consists of determination of 
the overall characteristics of the compressor, velocity distributions at 
the inlet of IGV and after each blade row, check for axisymmetry, check 
for unsteadiness at the inlet. Wherever possible, these results are com- 
pared with those of Smith (13). In the following sections, a brief 
summary of the results and interpretation is given. 

Overall characteristics: The compressor was tested extensively with 

various configurations both with the auxiliary fan and without auxiliary 
fan, to get the maximum pressure rise. It was decided to retain the 
auxiliary fan permanently in the facility. All tests were carried out 
with the auxiliary fan operating at constant speed and the compressor 
rotor operating at 1100 r.p.m. The throttle setting was varied to deter- 
mine the performance characteristics of the compressor at various flow 
coefficiants. The procedure for these tests was as follows. 

For each throttle setting, total pressure, static pressure and flow 
angle traverses were taken using pitot probe and two-hole wedge probe, at 
four axial stations 1, 2, 3 and 5. The wall static pressures were also 
taken at the hub and the casing to check the validity of static pressure 
traverses. Radial distribution of axial and tangeptial velocities and 
stagnation pressure coefficient were determined from the measured data. 

By integrating over the annulus, the mass-averaged pressure-rise coeffi- 
cient ^ was determined. The flow coefficient $ was determined by integrat- 
ing the axial velocity profile at the IGV inlet over the annulus area. 

Similarly $ was determined at other arial stations and compared with 
6 at the inlet of IGV. For the four runs that have been made, the flow 

*These and other notations are defined at the end of this section. 
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coefficients after 1GV, rotor and stator agree within 4Z with the flow 
coefficient determined from velocity profile at IGV inlet. 

A plot of the pressure-rise coefficient ^ vs. flow coefficient $ is 
shown in Fig. 5. With the present arrangement, the throttle cannot be 
opened any further than that corresponding to the highest capacity data 
point in the figure. Smith's data was also plotted in the same figure. 

In general, the present curve is lower than that of Smith by about 10Z. 

Radial distribution of flow properties: In Fig. 6, the tangential 
and axial velocity distributions, flow angle and local pressure-rise coef- 
ficient distributions for $ - 0.516 (design) are shown. Also shown in the 
same figure are the design values as well as the measured values of Smith 
(13). Agreement between the radial distribution of axial velocity in 
both cases is good, the major discrepancy being observed at the hub and 
tip. The tangential velocity and flow angle distributions are about 10Z 
lower than those measured by Smith and 5Z lower than design values. 

In the same figure, design distribution of axial and tangential 
velocities are shown. These distributions are obtained by axlsynmetric 
design method which does not include secondary flow effects. Smith 
observed that the difference between his measured tangential velocity 
distribution and design distribution is due to the over-turning in the 
guide vanes. 

In Fig. 7 radial distributions of relative flow quantities at rotor 
exit at $ ■ 0.516 for the present configuration, measured quantities by 

w e w 

Smith (13) and design values are presented. In general, , rr- and B 

u t t 

for the present case are higher than those obtained by Smith. This is to 

v e 

be expected, since is lower for the present case. 
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Design and measured velocity diagrams at six radial locations for the 
compressor stage is shown in Fig. 8. For inlet guide vanes , the agree- 
ment between design and measured velocity diagrams is good except at the 
hub and tip. This may be due to the secondary flow, which is not included 
in the design velocity diagram. The agreement between design and measured 
velocity triangles at rotor inlet Is poor. The reason for this discrepancy 
is not evident. The agreement between design and measured velocity 
diagrams at rotor is good except at the tip. Here the difference is con- 
siderable. This may be due to tip leakage flows. 

Stator inlet angles are those measured at the rotor exit. Hence the 
design and measured velocity diagrams are in good agreement. The spacing 
between rotor and stator is more than that of Smith (13). Hence the stator 
inlet flow is different from that of Smith. So the stator exit velocity 
diagrams are different from the design velocity diagrams. 

Axis ymme try and unsteadiness test: Results of this test are shown in 

Fig. 9. Radial surveys of the mean velocity and turbulence intensities 
were measured by a single hot-wire at inlet to the 1GV at three circum- 
ferential stations, shown in the insert of Fig. 9. 

Turbulence intensities at stations 1 and 2 were in good agreement and 
usually less than 2%. Turbulence intensities at station 3 are higher than 
those at the other two stations. This may be due to the proximity of 


station 3 to the floor 
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Noaenclature for Chap tar IV 
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radius 
huh radius 
rotor tip radius 

total pressure rise across the rotor 

rotor tip velocity 

relative velocity 

absolute flow angle 

relative flow angle 

mass averaged flow coefficient 



local pressure rise coefficient 



mass averaged pressure rise coefficient 
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Subscripts 

1 inlet to a blade row 

2 exit of a blade row 

x axial direction 


e 


tangential direction 
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Fig. 7. Radial Distribution of Wq/Uj., W/U t and f? 

at Rotor Exit at $ - 0.516 and x ■ 0.466. 
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